Nesting of electron and hole Fermi surfaces in non-superconducting BaFe2P2 
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We report detailed measurements of the de Haas- van Alphen (dHvA) effect in BaFe2P2 which is 
the end member of the superconducting series, BaFe2(Asi_a;Pi;)2. Using high purity samples we 
are able to observe dHvA oscillations from all the sheets of Fermi surface and hence build up a 
complete detailed picture of its structure. The results show the existence of a highly warped section 
of hole surface which may be the origin of the nodal gap structure in this series. Importantly, we 
find that, even for this non-superconducting end member, one of the hole surfaces almost exactly 
nests with the inner electron surface. This suggests that improved geometric nesting does not drive 
the increase in pairing strength with decreasing x in this system. 



o 
o 

Oh 
=3 

^ ! 

S; 

I 

x3 : 

a • 
o ■ 
o 



> 
o 

o 

^' 

o 



X 



A common feature that links together the various fam- 
ilies of iron-based superconductors is their unusual elec- 
tronic structure, vifhich consists of almost nested, quasi 
two dimensional electron and hole Fermi surfaces (FS). 
This quasi nesting produces strong peaks in the real and 
imaginary parts of the susceptibility (both bare and in- 
teracting), which in many theories [K, Q drives the su- 
perconductivity via a spin-fluctuation mechanism. How- 
ever, the relative importance of geometric nesting and the 
many-body interactions {U and J in Hubbard type mod- 
els) remains under debate. Studies of how the FS evolves 
and correlates with the changing material properties over 
the phase diagram are therefore of great importance. 

Although, quasi nested electron and hole bands are a 
common feature of calculated band-structures of these 
compounds, the exact topology and size of the FS is 
very difficult to compute from first principles. There 
are usually four or five bands which cross the Fermi en- 
ergy, giving rise to small Fermi pockets each of which fills 
~ 10 — 20% of the Brillouin zone volume. The small size 
of these pockets means their size and topology are very 
sensitive to small shifts in the band energies, which can 
result from small structural changes or from interactions 
with spin-fluctuations which scatter electrons between 
the electron and hole sheets 0]. Hence, it is important 
to experimentally determine the full three dimensional 
FS so that quantitative tests of the theory can be made. 
Measurements of the de Haas-van Alphen (dHvA) effect 
provide a very accurate way of measuring the bulk FS. 

A particularly useful system for studying the influence 
of FS topology on the ground state of the iron-based su- 
perconductors is BaFe2(Asi_a;P2;)2. Here, the Neel tem- 
perature of spin density wave ground state of the end 
member BaFe2As2 is driven to zero and superconductiv- 
ity emerges with a maximum Tc = 30 K, as As is sub- 
stituted with the isovalent element P. So unlike charge 
doped iron-pnictide systems, (e.g., Bai_a;Ka;Fe2As2), 
BaFe2(Asi_3;Pj;)2 is expected to remain perfectly com- 
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FIG. 1. (a) An example torque versus field trace with 
background subtracted (b) Fast Fourier transforms of the 
t{H) data (from 6-18 T)at various angles showing the dis- 
tinct dHvA peaks. Harmonics are identified with the prefix 
2. All data was taken at base temperature T ~ 0.35 K. Note 
the peaks below 0.5 kT are due to noise and features in the 
background cantilever magnetoresistance. 



pensated (equal volumes of electron and hole pockets) 
for all X. Minimal disorder is introduced in the Fe plane 
by these substitutions so that dHvA oscillations can be 
observed over a wide range of 2; (0.42 < a; < 1) [4]. 

Structurally, increasing x results in a decrease in the 
c/a ratio and a movement of the pnictogen atom towards 
the Fe plane. Density functional theory (DFT) calcula- 
tions based on the local density approximation (and its 
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FIG. 2. (color online) dHvA frequencies multiplied by cosine of the rotation angle versus angle, (a) Shows the calculations for 
the unshifted bands rotating from [001] towards [100]. In (b) the symbols show experimental data for rotations from [001] to 
both [100] and [110] whereas the solid lines show the calculations for the energy shifted bands as described in the text (note 
that the predicted orbit Si which has F ~ 100 T at = 0° is not shown) (c) Shows the origin of the orbits on the calculated 
electron (e) and hole (h) FS sheets at an angle of 9 = 20° (towards [110]). 



variants) predict that throughout the series the electron 
sheets change very little, but for the hole sheets there are 
much greater changes [4]. BaFe2As2 has three quasi 2D 
hole sheets but for BaFe2P2 one of these disappears and 
one of the hole sheets develops pronounced warping close 
to the top of the zone (dumbbell, band 2 in Fig. A 
notable property of the BaFe2(Asi_j:Pa;)2 series is that 
there is clear evidence from multiple probes [1, 0] that 
there are line nodes in its superconducting energy gap. 
It has been proposed theoretically that these nodes are 
located on the dumbbell hole surface [s'] although there 
is experimental evidence that rather points to the nodes 
being located on the electron sheets 0, . Establishing 
a clear connection between the unusual gap structure of 
BaFe2(Asi_xPa;)2 and its FS topology would be a strong 
confirmation that the pairing mechanism is indeed medi- 
ated by spin-fluctuations in these materials. 

In previous studies, dHvA signals from all of the FS 
sheets were not observed and so the complete topology 
could not be determined. In particular, Shishido et al. [4] 
only observed signals from the electron sheets. Follow- 
ing this, Analytis et al. fll| reported measurements for 
X — 0.63 where orbits on the smaller hole sheet were ob- 
served and were found to correspond to almost the same 
crosssectional area as the smaller electron pocket. 

Here we report new measurements of the dHvA effect 
for the end member BaFe2P2. The much higher purity 
level of our new samples allows us to observe dHvA os- 
cillations originating from all the FS sheets and hence 
we are able to make a precise determination of the com- 
plete FS. A surprising aspect of the results is that they 
show that even for this non-superconducting (and non- 
magnetically ordered) composition, one of the hole sheets 
is almost perfectly nested with the (smaller) electron FS, 



similar to that found for the superconducting composi- 
tion X = 0.63 This suggests that nesting is probably 
a necessary but not sufficient ingredient for superconduc- 
tivity. 

Samples of BaFe2P2 were grown as described in Ref. 
[l^ . The lattice parameters and the crystal orientation 
was determined by x-ray diffraction. The best results 
were obtained on a sample with dimensions 50 x 75 x 
20/im, which was mounted on micro piezo resistive can- 
tilever in order to measure the magnetic torque. The 
torque sensor was fixed to a rotating probe, inside a "^He 
cryostat equipped with a superconducting magnet. The 
orientation of the sample on the cantilever was deter- 
mined by comparing optical images of the mounted sam- 
ple with images taken on the x-ray mount. We estimate 
the accuracy of the in-plane alignment 4> is around ±5°. 

Example torque r versus field H data is shown in 
Fig. [T] In this sample, the oscillations are visible down 
to fields as low as 5T at some angles. By taking fast 
Fourier transforms of the t{B^^) data we identified 8 
different dHvA orbits. As many of these peaks are very 
close in frequency they can only be clearly identified by 
performing long field sweeps (typically between 6 and 
18 T) with small (1 degree) angular increments. In this 
way the various crossing branches can be easily identi- 
fied. Complete rotations were performed in two principle 
directions, namely from H perpendicular to the Fe plane 
[001] towards either the [100] or [110] directions. In both 
cases we define the rotation angle 6 ~ {) for i7||[001]. 

For a perfectly two dimensional FS sheet, the dHvA 
frequency F (x 1/ cos{9), hence for quasi- two dimensional 
systems the dHvA data are best appreciated by plotting 
Fcos6 versus 9 as in Fig. [5] In such a plot, a local 
minimum in the FS crosssection (neck) gives rise to a 
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FIG. 3. (color online) Crosssections of the determined Fermi surface of BaFe2P2 in: (a) the ab plane at three different fez values 
(quoted in units of c*), (b) a c-axis cut along the [110] direction and through the F point. The largest band character element 
at each fc-point is indicated [13 ]. 



branch where F cos Q increases with 6*, whereas for a local 
maximum (belly) F cos 6 decreases with d. 

In order to make a more precise analysis of the re- 
sults we compare the data with predicted extremal orbits 
derived from bandstructure calculations. These density 
functional calculations were done using an augmented 
plane wave plus local orbits scheme as implemented in 
the Wien2K package [3, 15 1. The calculated extremal 
orbits as function of the field angle for a [001] to [100] ro- 
tation are shown in Fig. [SJa) . Although this calculation 
reproduces the general features of the experiment, the 
exact topology and FS size is not well reproduced. As 
we have found in other pnictide systems 3 12 1 the cal- 
culation can be brought into much better agreement with 
experiment by performing rigid (fc-independent) shifts of 
the band energies. The two electron bands (4 and 3) 
are easily identified as the a and /3 branches respectively 
and these are brought into almost exact agreement by 
shifting the energies of band 4 up by 68 meV and band 
3 up by 58 meV. These surfaces have quite a compli- 
cated topology which is well reproduced in the calcula- 
tions [see Fig. Ellb)]. Note the large differences between 
the rotation toward [100] and [110]. The large (high- 
est fundamental frequency) orbit 72 originates from the 
highly warped part of the band 2 hole 'dumbbell' surface 
centered around the top of the zone. The calculations 
reproduce this branch without any applied energy shift. 
The curvature of Fcos(0) for the 81 branch clearly iden- 
tifies it as a maximum (belly) crosssection and so it must 
originate from band 1 as this is the only remain uniden- 
tified maximal orbit. Shifting the energy of this band 
down by 113 meV gives an almost exact match. With 
this shift, the minimal crosssection Si of this band would 
give a dHvA frequency around 100 T. We do not see con- 
clusive evidence for this orbit probably because its ex- 
pected weak signal is likely buried in the low frequency 
cantilever noise (see Fig.[T]). Finally, we can identify the 
last branch 71, as the minimum of band 2. Although 
no energy shift was required to explain the maximum of 
this band, we require a shift down of 52 meV to fit the 
minimum. With these energy shifts, the FS is accurately 



TABLE I. Measured effective masses m* (uncertainty ~2%) 
along with mass enhancements m* /mt calculated using the 
bandstructure values (rrii,) from the energy shifted bands. The 
mean free paths £ are also reported [l^ . 
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determined. As a check we calculate the number of elec- 
tron and holes in each shifted FS sheet and find that, 
as expected, they are equal to within 0.017 electrons per 
unit cell. 

It has been suggested that physically these shifts result 
from a spin fluctuation interaction [3] . Scattering of elec- 
trons by the spin-fluctuations between the electron and 
hole surfaces causes both surfaces to shrink in volume 
and hence the energy shifts are in opposite directions. 
This is supported by results for CaFe2P2 ^18^, which has 
a non-nested FS similar to the collapsed tetragonal, high 
pressure phase of CaFe2As2, and where the DFT calcu- 
lations are in excellent agreement with experiment with- 
out any shifting. Further, for BaFe2(Asi_a;Pa;)2 it was 
found that the size of energy shifts varied approximately 
linearly with x, becoming larger as the SDW phase tran- 
sition is approached [4,] and the spin-fluctuations grow in 
magnitude. Our results here also support this interpre- 
tation, because we find that most sections of the FS of 
BaFe2P2 are close to being geometrically nested and are 
smaller compared to the DFT calculation whereas the 
highly warped part of band 2 dumbbell (72 orbit) is far 
from nesting and remains unchanged. This supports the 
suggestion of Suzuki et al. Q that this part of the Fermi 
surface does not contribute strongly to the pairing and is 
therefore may be susceptible to gap node formation. 

The masses m* of the quasiparticles on each of the 
extremal orbits were determined by measuring the tem- 
perature dependence of the dHvA oscillation amplitude 
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|19| . Fits to this data to determine m* are reported in 
Table. HI along with the masses mt calculated from the 
shifted band-structure bands for the same orbits and at 
the same field angles. The mass enhancements m*/mi, 
are all relatively similar and small ranging from 1.6 to 
1.9. Note the values for the /3 orbits are smaller than the 
value give in Ref. Q but consistent within the error. 

An interesting feature of the results is that they show 
that for this non-superconducting end member of the 
BaFe2(Asi_^Pj;)2 series, one of the electron and hole 
sheets are very well nested. The closeness of the crosssec- 
tional areas of the inner electron and hole FS is evident 
from the rotation plots shown in Fig. [2] The fact that 
these areas are so close explains why they were not ob- 
served in previous studies, where the dHvA oscillations 
were not observable over such a wide range of inverse field 
as that were here. In Fig.|3]we show sections through the 
determined FS. Although the in-plane area of the elec- 
tron sheet FS does not change strongly with , its shape 
does. Because of the symmetry of the body centered 
tetragonal unit cell the crosssection at the zone center 
is exactly the same as at the zone top but rotated by 
7r/2. The best match, where the inner electron sheet is 
almost perfectly circular is half way between these points 
at kz = 0.25. Furthermore the band character of these 
two sheets is also very similar; they both have predomi- 
nately Fe d^z/yz character. 

Comparing our data with that of Analytis et al. 11] for 
the superconducting (Tc = 7K) composition x — 0.63 
suggests that besides the above mentioned shrinking of 
the overall size of the Fermi surface there is remark- 
ably little change in the overall topology with x (for 
0.63 < X < 1). A close comparison between the two 
data sets [3| shows that the F{9) data may be almost 
perfectly superimposed by simply scaling the mean fre- 
quencies of the electron bands. This suggests that the 
c-axis warping of the electron and hole sheets remain re- 
markably constant as a function of P content x. The 
anomalous downturn in the frequency of the electron a 
branch between 6 = 40° and 50° identified in Ref. 11 1 
was correctly assigned to a crossing of the hole and elec- 
tron orbit frequencies. However, the higher scattering in 
these doped samples probably prevented the observation 
of the dumbbell (72) orbit. In Ref. ll| the 7 branch was 
assigned to the minimum of band 1. However, as the 
size and dispersion of this branch is almost identical to 
the 71 branch in the present work we suggest that this 
assignment was incorrect, and it is instead the minimum 
of band 2 as in the present work. This effects somewhat 
the fitted band structure of Ref. [llj but does not change 
the conclusion that the electron and hole bands are well 
nested for x = 0.63. 

Combined with the present results it appears then that 
rather than the bands become much better nested as x 
approaches that optimal for superconductivity, the rela- 
tive size of the electron and hole sheets remain remark- 



ably constant. As the electron sheets shrink with decreas- 
ing x the hole bands shrink by a similar amount with- 
out any major change in topology. This is not expected 
from conventional bandstructure calculations, where the 
hole band topology changes very dramatically as x is de- 
creased. It is therefore likely that the dumbbell shape 
of the larger hole band remains for the superconduct- 
ing compositions. Recent photoemission results for the 
X — 0.38 composition [i^l indicate the FS topology of 
this close to optimal Tc composition is indeed close to 
that observed here for the end member. 

In summary, our results determine, in very high de- 
tail, the structure of the FS of the end member of 
the BaFe2(Asi_a;P2;)2 series. Together with the known 
shrinking of the average size of the electron pockets and 
invoking the expected compensated nature of the FS, it 
is now possible to estimate the topology of both the elec- 
tron and hole surface across the phase diagram. An im- 
portant aspect of the results is that they show that one of 
the hole and electron FS are well nested even for this non- 
superconducting (and non magnetic) composition. This 
suggests that the increase in pairing susceptibility as x 
is decreased is caused by increases in the many body in- 
teraction parameters U and J rather than any effect of 
improved geometric nesting. 
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